    local controls lnsize age1900 fxassratio_MEAN
    local depvar DR
    local varlist1 Ataxpayers_s `controls'
    local varlist2 Acivil_s `controls'
    local varlist3 Acivil_s Awealth_s `controls'
    local varlist4 Acivil_s Awealth_s Doslo `controls'

    logit `depvar' `varlist1' , nolog robust
    margins, at(Ataxpayers_s=(.2587 .2874)) asobserved post             /* one std dev around the mean */
    di   _b[1._at] - _b[2._at]
    test _b[1._at] = _b[2._at]

    qui logit `depvar' `varlist1' , nolog robust
    margins, at(lnsize=(.754 1.377)) asobserved post             /* one std dev around the mean */
    di   _b[1._at] - _b[2._at]
    test _b[1._at] = _b[2._at]

    qui logit `depvar' `varlist1' , nolog robust
    margins, at(age1900=(37.43 64.26)) asobserved post             /* one std dev around the mean */
    di   _b[1._at] - _b[2._at]
    test _b[1._at] = _b[2._at]

    qui logit `depvar' `varlist1' , nolog robust
    margins, at(fxassratio_MEAN=(.312 .488))  asobserved post             /* one std dev around the mean */
    di   _b[1._at] - _b[2._at]
    test _b[1._at] = _b[2._at]

    logit `depvar' `varlist2' , nolog robust
    margins, at(Acivil_s=(.0152   .0217)) asobserved post             /* one std dev around the mean */
    di   _b[1._at] - _b[2._at]
    test _b[1._at] = _b[2._at]

    qui logit `depvar' `varlist2' , nolog robust
    margins, at(lnsize=(.754 1.377)) asobserved post             /* one std dev around the mean */
    di   _b[1._at] - _b[2._at]
    test _b[1._at] = _b[2._at]

    qui logit `depvar' `varlist2' , nolog robust
    margins, at(age1900=(37.43 64.26)) asobserved post             /* one std dev around the mean */
    di   _b[1._at] - _b[2._at]
    test _b[1._at] = _b[2._at]

    qui logit `depvar' `varlist2' , nolog robust
    margins, at(fxassratio_MEAN=(.312 .488))  asobserved post             /* one std dev around the mean */
    di   _b[1._at] - _b[2._at]
    test _b[1._at] = _b[2._at]


    logit `depvar' `varlist3' , nolog robust
    margins, at(Acivil_s=(.0152   .0217)) asobserved post             /* one std dev around the mean */
    di   _b[1._at] - _b[2._at]
    test _b[1._at] = _b[2._at]

    qui logit `depvar' `varlist3' , nolog robust
    margins, at(Awealth_s=(1.225 1.680))  asobserved post             /* one std dev around the mean */
    di   _b[1._at] - _b[2._at]
    test _b[1._at] = _b[2._at]

    qui logit `depvar' `varlist3' , nolog robust
    margins, at(lnsize=(.754 1.377)) asobserved post             /* one std dev around the mean */
    di   _b[1._at] - _b[2._at]
    test _b[1._at] = _b[2._at]

    qui logit `depvar' `varlist3' , nolog robust
    margins, at(age1900=(37.43 64.26)) asobserved post             /* one std dev around the mean */
    di   _b[1._at] - _b[2._at]
    test _b[1._at] = _b[2._at]

    qui logit `depvar' `varlist3' , nolog robust
    margins, at(fxassratio_MEAN=(.312 .488))  asobserved post             /* one std dev around the mean */
    di   _b[1._at] - _b[2._at]
    test _b[1._at] = _b[2._at]


    logit `depvar' `varlist4' , nolog robust
    margins, dydx(`varlist4')  asobserved

    qui logit `depvar' `varlist4' , nolog robust
    margins, at(Acivil_s=(0.0153   0.0217)) asobserved post             /* one std dev around the mean */
    di   _b[1._at] - _b[2._at]
    test _b[1._at] = _b[2._at]

    qui logit `depvar' `varlist4' , nolog robust
    margins, at(Awealth_s=(1.2481 1.6882))  asobserved post             /* one std dev around the mean */
    di   _b[1._at] - _b[2._at]
    test _b[1._at] = _b[2._at]

    qui logit `depvar' `varlist4' , nolog robust
    margins, at(lnsize=(.754 1.377)) asobserved post             /* one std dev around the mean */
    di   _b[1._at] - _b[2._at]
    test _b[1._at] = _b[2._at]

    qui logit `depvar' `varlist4' , nolog robust
    margins, at(age1900=(37.43 64.26)) asobserved post             /* one std dev around the mean */
    di   _b[1._at] - _b[2._at]
    test _b[1._at] = _b[2._at]

    qui logit `depvar' `varlist4' , nolog robust
    margins, at(fxassratio_MEAN=(.312 .488))  asobserved post             /* one std dev around the mean */
    di   _b[1._at] - _b[2._at]
    test _b[1._at] = _b[2._at]

